Magnetic catalysis and anisotropic confinement in QCD 
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The expressions for dynamical masses of quarks in the chiral limit in QCD in a strong magnetic 
field are obtained. A low energy effective action for the corresponding Nambu-Goldstone bosons 
is derived and the values of their decay constants as well as the velocities are calculated. The 
existence of a threshold value of the number of colors A^*'"^, dividing the theories with essentially 
different dynamics, is established. For the number of colors Nc <C A^*'"^, an anisotropic dynamics of 
confinement with the confinement scale much less than Aqcd and a rich spectrum of light glueballs 
is realized. For Nc of order A^*'"' or larger, a conventional confinement dynamics takes place. It 
is found that the threshold value Af^'*' grows rapidly with the mag netic field [A^^'"' > 100 for 
\eB\ > (1 GeV)^]. In contrast to QCD with a nonzero baryon density, there are no principal 
obstacles for examining these results and predictions in lattice computer simulations. 
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I. INTRODUCTION 

Since the dynamics of QCD is extremely rich and complicated, it is important to study this theory under external 
conditions which provide a controllable dynamics. On the one hand, this allows one to understand better the vacuum 
structure and Green's functions of QCD, and, on the other hand, there can exist interesting applications of such 
models in themselves. The well known examples are hot QCD (for a review see Ref. [Q) and QCD with a large baryon 
density (for a review see Ref. [||). 

Studies of QCD in external electromagnetic fields had started long ago [^,||. A particularly interesting case is an 
external magnetic field. Using the Nambu-Jona-Lasinio (NJL) model as a low energy effective theory for QCD, it 
was shown that a magnetic field enhances the spontaneous chiral symmetry breakdown. The understanding of this 
phenomenon had remained obscure until a universal role of a magnetic field as a catalyst of chiral symmetry breaking 
was established in Refs. [||J^. The general result states that a constant magnetic field leads to the generation of a 
fermion dynamical mass (i.e., a gap in the one-particle energy spectrum) even at the weakest attractive interaction 
between fermions. For this reason, this phenomenon was called the magnetic catalysis. The essence of the effect is 
the dimensional reduction D —^ D — 2 in the dynamics of fermion pairing in a magnetic field. In the particular case of 
weak coupling, this dynamics is dominated by the lowest Landau level (LLL) which is essentially D — 2 dimensional 
1^,^ . The applications of this effect have been considered both in condensed matter physics ||^,^ and cosmology (for 
reviews see ReL [^). 

The phenomenon of the magnetic catalysis was studied in gauge theories, in particular, in QED ||ic|-[l5[| and in QCD 
p6|-p^. In the recent work it has been suggested that the dynamics underlying the magnetic catalysis in QCD is 
weakly coupled at sufficiently large magnetic fields. In this paper, we study this dynamical problem rigorously, from 
first principles. In fact, we show that, at sufficiently strong magnetic fields, \eB\ ^ ^qcd' there exists a consistent 
truncation of the Schwinger-Dyson (gap) equation which leads to a reliable asymptotic expression for the quark mass 



niq. Its explicit form reads: 
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where is the electric charge of the q-th quark and Nc is the number of colors. The numerical factors Ci and C2 
equal 1 in the leading approximation that we use. Their value, however, can change beyond this approximation and 
we can only say that they are of order 1. The constant c„ is defined as follows: 
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where iV„ and Nd are the numbers of up and down quark flavors, respectively. The total number of quark flavors is 
Nf = Nu + Nd- The strong coupling as in the last equation is related to the scale y^eSf, i.e., 

— ~ b\n ' ' , where b= i. (3) 

^QCD 



We should note that in the leading approximation the energy scale •y/|ei?| in Eq. (g) is fixed only up to a factor of 
order 1. 

As we discuss below, because of the running of a^, the value of the dynamical mass grows very slowly with 
increasing the value of the background magnetic field. Moreover, there may exist an intermediate region of fields 
where the mass decreases with increasing the magnetic field. Another, rather unexpected, consequence is that a strong 
external magnetic field suppresses the chiral vacuum fluctuations leading to the generation of the usual dynamical 
mass of quarks rnf^^ ~ 300 MeV in QCD without a magnetic field. In fact, in a wide range of strong magnetic fields 
A'^ ^ -B ^ (10 TeV)^ (where A is the characteristic gap in QCD without the magnetic field; it can be estimated to be 



a few times larger than Aqcd), the dynamical mass dll) remains smaller than rrSdvn- ^^^^ shown in Sec. IV, 



this point is intimately connected with another one: in a strong magnetic field, the confinement scale, Xqcd, is much 
less than the confinement scale Aqcd in QCD without a magnetic field. 

The central dynamical issue underlying this dynamics is the effect of screening of the gluon interactions in a magnetic 
field in the region of momenta relevant for the chiral symmetry breaking dynamics, |fc^| ^ In this region, 

gluons acquire a mass Mg of order Nfas\eqB\. This allows to separate the dynamics of the magnetic catalysis from 
that of confinement. More rigorously, Mg is the mass of a quark-antiquark composite state coupled to the gluon field. 
The appearance of such mass resembles pseudo-Higgs effect in the 1 + 1 dimensional massive QED (massive Schwinger 
model) ||l^ (see below). 

Since the background magnetic field breaks explicitly the global chiral symmetry that interchanges the up and down 
quark flavors, the chiral symmetry in this problem is SU{Nu)l x SU{nJ)r x SU{Nd)L x SU{Nd)R x U'^-^1)a- The 
L/'^^'(l)yi is connected with the current which is an anomaly free linear combination of the U'^'^\1)a and [/'^"^(l)yi 
currents. [The U^~\1)a symmetry is of course absent if either Nd or N^ is equal to zero]. The generation of quark 
masses breaks this symmetry spontaneously down to SU{Nu)v x SU{Nd)v and, as a result, N^ + N^ — 1 gapless 



Nambu-Goldstonc (NG) bosons occur. In Sec. Ill, we derive the effective action for the NG bosons and calculate 
their decay constants and velocities. 

The present analysis is heavily based on the analysis of the magnetic catalysis in QED done by Gusynin, Miransky, 
and Shovkovy A crucial difference is of course the property of asymptotic freedom and confinement in QCD. In 
connection with that, our second major result is the derivation of the low energy effective action for gluons in QCD 
in a strong magnetic field [see Eq. (|l8|) below]. The characteristic feature of this action is its anisotropic dynamics. 
In particular, the strength of static (Coulomb like) forces along the direction parallel to the magnetic field is much 
larger than that in the transverse directions. Also, the confinement scale in this theory is much less than that in QCD 
without a magnetic field. This features imply a rich and unusual spectrum of light glueballs in this theory. 

A special and interesting case is QCD with a large number of colors, in particular, with TV, ^ oo (the 't Hooft limit). 
In this limit, the mass of gluons goes to zero and the expression for the quark mass becomes essentially different [see 
Eq. ( p5| ) in Sec. Q. In fact, it will be shown that, for any value of an external magnetic field, there exists a threshold 
value A^*''^ rapidly growing with \eB\ [e.g., A^*'"^ > 100 for \eB\ > (1 GeV)^]. For Nc of the order N*'''' or larger, the 
gluo n mass becomes small and irrelevant for the dynamics of the generation of a quark mass. As a result, expression 
( psi ) for ruq takes place for such large Nc- The confinement scale in this case is close to Aqcd- Still, as is shown in 
Sec. ^ the dynamics of chiral symmetry breaking is under control in this limit if the magnetic field is sufficiently 
strong. 

It is important that, unlike the case of QCD with a nonzero baryon density, there are no principal obstacles for 
checking all these results and predictions in lattice computer simulations of QCD in a magnetic field. 

II. MAGNETIC CATALYSIS IN QCD 

We begin by considering the Schwinger-Dyson (gap) equation for the quark propagator. It has the following form: 
G-\x, y) = S-^{x, y) + 47ra,7'' / G{x, z)T''{z, y, z')V,^{z', x)d^zd^z', (4) 



where S{x, y) and G{x^ y) are the bare and full fermion propagators in an external magnetic field, ^^^{x, y) is the full 
gluon propagator and V^ix, y, z) is the full amputated vertex function. Since the coupling as related to the scale |ei3| 
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is small, one might think that the rainbow (ladder) approximation is reliable in this problem. However, this is not the 
case. Because of the (l+l)-dimensional form of the fermion propagator in the LLL approximation, there are relevant 
higher order contributions . Fortunately one can solve this problem. First of all, an important feature of the 

quark-antiquark pairing dynamics in QCD in a strong magnetic field is that this dynamics is essentially abelian. This 
feature is provided by the form of the polarization operator of gluons in this theory. The point is that the dynamics 
of the quark-antiquark pairing is mainly induced in the region of momenta k much less than •\/|ei3| . This implies that 
the magnetic field yields a dynamical ultraviolet cutoff in this problem. On the other hand, while the contribution of 
(electrically neutral) gluons and ghosts in the polarization operator is proportional to fc^, the fermion contribution is 
proportional to \eqB\ As a result, the fermion contribution dominates in the relevant region with <C |ei?|. 

This observation implies that there are three, dynamically very different, scale regions in this problem. The first 
one is the region with the energy scale above the magnetic scale ^y\eB\. In that region, the dynamics is essentially 
the same as in QCD without a magnetic field. In particular, the running coupling decreases logarithmically with 
increasing the energy scale there. The second region is that with the energy scale below the magnetic scale but much 
larger than the dynamical mass m^. In this region, the dynamics is abelian like and, therefore, the dynamics of the 
magnetic catalysis is similar to that in QED in a magnetic field. At last, the third region is the region with the energy 
scale less than the gap. In this region, quarks decouple and a confinement dynamics for gluons is realized. 

Let us first consider the intermediate region relevant for the magnetic catalysis. As was indicated above, the 
important ingredient of this dynamics is a large contribution of fermions to the polarization operator. It is large 
because of an (essentially) 1+1 dimensional form of the fermion propagator in a strong magnetic field. Its explicit 
form can be obtained by modifying appropriately the expression for the polarization operator in QED in a magnetic 
field (l|: 

Nf 
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where g'^'^ = diag(l, 0, 0, — 1) is the projector onto the longitudinal subspace, and fc||' = 5||"'fcy (the magnetic field 
is in the direction). Similarly, we introduce the orthogonal projector = g'^" — g'^^'^ = diag(0, — 1, — 1, 0) and 

fc^ = gj^ki^ that we shall use below. Notice that quarks in a strong magnetic field do not couple to the transverse 
subspace spanned by g'^ and fc^ . This is because in a strong magnetic field only the quark from the LLL matter and 
they couple only to the longitudinal components of the gluon field. The latter property follows from the fact that 
spins of the LLL quarks arepolarized along the magnetic field pO[ |. 

The expressions (^ and (j6) coincide with those for the polarization operator in the massive Schwinger model if the 
parameter as\eqB\/2 here is replaced by the dimensional coupling ai of QEDi^i. As in the Schwinger model, Eq. (|^) 
implies that there is a massive resonance in the — ^f\9\\'^ component of the gluon propagator. Its mass is 

Nf 
q=l 

This is reminiscent of the pseudo-Higgs effect in the (l+l)-dimensional massive QED. It is not the genuine Higgs 
effect because there is no complete screening of the color charge in the infrared region with |fc^| ^ m^. This can be 
seen clearly from Eq. (^. Nevertheless, the pseudo-Higgs effect is manifested in creating a massive resonance and 
this resonance provides the dominant forces leading to chiral symmetry breaking. 

Now, after the abelian like structure of the dynamics in this problem is established, we can use the results of the 
analysis in QED in a magnetic field |pT[| by introducing appropriate modifications. The main points of the analysis 
are: (i) the so called improved rainbow approximation is reliable in this problem provided a special non-local gauge 
is used in the analysis, and (ii) for a small coupling Us (a in QED), the relevant region of momenta in this problem is 
<^ <C \eB\. We recall that in the improved rainbow approximation the vertex T'^{x,y,z) is taken to be bare 
and the gluon propagator is taken in the one-loop approximation. Moreover, as we argued above, in this intermediate 
region of momenta, only the contribution of quarks to the gluon polarization tensor (|6|) matters. [It is appropriate to 
call this approximation the "strong-magnetic- field- loop improved rainbow approximation" . It is an analog of the hard- 
dense-loop improved rainbow approximation in QCD with a nonzero baryon density |po| ]. As to the modifications, 
they are purely kinematic: the overall coupling constant in the gap equation a and the dimensionless combination 
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M'^/\eB\ in QED have to be replaced by asiN"^ - l)/2Nc and M^/\eqB\, respectively. This leads us to the expression 
(|l|) for the dynamical gap. 

After expressing the magnetic field in terms of the running coupling, the result for the dynamical mass takes the 
following convenient form: 
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a,(iV2_i)in(C2/c,a,) 



(8) 



As is easy to check, the dynamical mass of the w-quark is considerably larger than that of the d-quark. It is also 
noticeable that the values of the w-quark dynamical mass becomes comparable to the vacuum value m|;°J„ ~ 300 MeV 
only when the coupling constant gets as small as 0.05. 

Now, by trading the coupling constant for the magnetic field scale \eB\, we get the dependence of the dynamical 
mass on the value of the external field. The numerical results are presented in Fig. [we used Ci = C2 = 1 in Eq. 
(I)]. 
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FIG. 1. The dynamical masses of quarks as functions of hi{\eB\/AQCj}) for Nc — 3 and two different values of Nf = Nu + Nd'. 
(i) masses of u-quark (solid line) and d-quark (dash-dotted line) for Nu = 1 and Nd = 2; (ii) masses of u-quark (dashed line) 
and d-quark (dotted line) for A''^ — 2 and Nd ~ 2. The result may not be reliable in the weak magnetic field region (shaded) 
where some of the approximations break. The values of masses are given in units of Aqcd ~ 250 MeV. 



As one can see in Fig. ^ the value of the quark gap in a wide window of strong magnetic fields, Aq;-^^ ^ \eB\ < 

(10 TeV)^, remains smaller than the dynamical mass of quarks m^^y^ ~ 300 MeV in QCD without a magnetic field. 
In other words, the chiral condensate is partially suppressed for those values of a magnetic field. The explanation of 
this, rather unexpected, result is actually simple. The magnetic field leads to the mass Alg (0) for gluons. In a strong 
enough magnetic field, this mass becomes larger than the characteristic gap A in QCD without a magnetic field (A, 
playing the role of a gluon mass, can be estimated as a few times larger than Aqcd). This, along with the property of 
the asymptotic freedom (i.e., the fact that as decreases with increasing the magnetic field), leads to the suppression 
of the chiral condensate . 

This point also explains why our result for the gap is so different from that in the NJL model in a magnetic field 
1^. Recall that, in the NJL model, the gap logarithmically (i.e., much faster than in the present case) grows with a 
magnetic field. This is the related to the assumption that both the dimensional coupling constant G = g/A^ (with 
A playing a role similar to that of the gluon mass in QCD), as well as the scale A do not dependent on the value of 
the magnetic field. Therefore, in that model, in a strong enough magnetic field, the value of the chiral condensate is 
overestimated. 

The picture which emerges from this discussion is the following. For values of a magnetic field \eB\ ^ A^ the 
dynamics in QCD should be qualitatively similar to that in the NJL model. For strong values of the field, however, it 
is essentially different, as was described above. This in turn suggests that there should exist an intermediate region 
of fields where the dynamical masses of quarks decreases with increasing the background magnetic field. 



III. EFFECTIVE ACTION OF NG BOSONS 



The presence of the background magnetic field breaks explicitly the global chiral symmetry that interchanges the 
up and down quark flavors. This is related to the fact that the electric charges of the two sets of quarks are different. 
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However, the magnetic field does not break the global chiral symmetry of the action completely. In particular, in the 
model with the up quark flavors and the down quark flavors, the action is invariant under the chiral symmetry 
SU{Nu)l X SU{Nu)r X SU{Nd)L X SU{Nd)R x The C/(~'(1)a is connected with the current which is an 

anomaly free linear combinaton of the C/'^'*'(l)^ and U''^^{1)a currents. [The U^~'>{1)a symmetry is of course absent 
if either Nd or Nu is equal to zero] . 

The global chiral symmetry of the action is broken spontaneously down to the diagonal subgroup SU{Nu)v x 
SU{Nd)v when dynamical masses of quarks are generated. In agreement with the Goldstone theorem, this leads to 
the appearance of + Nj — 1 number of the NG gapless excitations in the low-energy spectrum of QCD in a strong 
magnetic field. Notice that there is also a pseudo-NG boson connected with the conventional (anomalous) U{1)a 
symmetry which can be rather light in a sufficiently strong magnetic field. 

Now, in the chiral limit, the general structure of the low energy action for the NG bosons could be easily established 
from the symmetry arguments alone. First of all, such an action should be invariant with respect to the space-time 
symmetry 5*0(1, 1) x 50(2) which is left unbroken by the background magnetic field [here the S0(1,1) and the S0(2) 
are connected with Lorentz boosts in the xq — 0:3 hyperplane and rotations in the xi ~ X2 plane, respectively]. Besides 
that, the low-energy action should respect the original chiral symmetry SU{Nu)l x SU {Nu)r x SU{Nd) l x SU{Nd) b. x 
C/^~^(l)^. These requirements lead to the following general form of the action: 

Cng ^ ^tr (^gi^^^d^E^d^El + i;^^^ 9^E„9,Et 
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(sf 9^Ea.Et + i'g^^^d^f:d,f:^) . (9) 



The unitary matrix fields = Gxp (^iJ2A=i^ ^^'^u / fu^ ^ ^rf = '^^P (J'J2a=i^ ■^^''^d / fdj ^ ^^'^ ^ = 6xp^iV27r// 

describe the NG bosons in the up, down, and U^^\1)a sectors of the original theory. The decay constants fu,fd,f 
and transverse velocities Vd, v can be calculated by using the standard field theory formalism (for a review, see for 
example the book [Q). Let us first consider the iV^ -|- — 2 NG bosons in the up and down sectors, assigned to the 
adjoint representation of the SU{Nu)v x SU{Nd)v symmetry. The basic relation is 



A 



(10) 



where = (^P'^,VgP±, P^^ and xf{k,P) is the Bethe-Salpeter (BS) wave function of the NG bosons {P is the 

momentum of their center of mass). In the weakly coupled dynamics at hand, one could use an analogue of the 
Pagels-Stokar approximation [p3| , p2| . In this approximation, the BS wave function is determined from the Ward 
identities for axial currents. In fact, the calculation of the decay constants and velocities of NG bosons resembles 
closely the calculation in the case of a color superconducting dense quark matter [p^ . In the LLL approximation, the 
final result in Euclidean space is 

The evaluation of this integral is straightforward. As a result, we get 



67r2 

Nr. 



fu-^\eBl (12) 



/^'^I^l^^l- (13) 

The remarkable fact is that the decay constants are nonzero even in the limit when the dynamical masses of quarks 
approach zero. The reason of that is the 1-1-1 dimensional character of this dynamics: as one can see from expression 
(|ll|), in the limit niq 0, the infrared singularity in the integral cancels the mass niq in the numerator. A similar 
situation takes place in color superconductivity [ p5| , p4[ : in that case the 1-1-1 dimensional character of the dynamics 
is provided by the Fermi surface. 

Notice that the transverse velocities of the NG bosons are equal to zero. This is also a consequence of the 1 -f- 1 
dimensional structure of the quark propagator in the LLL approximation. The point is that quarks can move in the 
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transverse directions only by hopping to higher Landau levels. Taking into account higher Landau levels would lead 



to nonzero velocities suppressed by powers of 



\eB\. In fact, the explicit form of the velocities was derived in 



the weakly coupled NJL model in an external magnetic field [see Eq. (65) in the second paper of Ref. |10)]. It is 



\mu,d\ 
\eB\ 



In- 



\eB\ 



\mu.d\ 



< 1. 



(14) 



A similar expression should take place also for the transverse velocities of the NG bosons in QCD. 

Now, let us turn to the NG boson connected with the spontaneous breakdown of the U^~\1)a- It is a SU{Nu)v x 
SU{Nd)v singlet. Neglecting the anomaly, we would actually get two NG singlets, connected with the up and down 
sectors, respectively. Their decay constants and velocities would be given by expression ( |l0| ) in which A"^ has to be 
replaced by A°. The latter is proportional to the imit matrix and normalized as the A"^ matrices: tr[(A°)^] = 2. It is 
clear that their decay constants and velocities would be the same as for the NG bosons from the adjoint representation. 
Taking now into account the anomaly, we find that the anomaly free U^~'^{1)a current is connected with the traceless 
matrix A°/2 
singlets. 



|1, d) and |1, w), as 



y/ Nu/NfXj)/2. Therefore the genuine NG singlet |1) is expressed through those two 

= y]Vd7]v7|i,w) 



y/Nu/Nf\l, d). This implies that its decay constant is 



[Ndfu + Nufdf {V2Nd + Nu)^N, 



\eB\ 



(15) 



Its transverse velocity is of course zero in the LLL approximation. 



IV. ANISOTROPIC CONFINEMENT OF GLUONS 

Let us now turn to the infrared region with |fc| ^ rud, where all quarks decouple (notice that we take here the 
smaller mass of d quarks). In that region, a pure gluodynamics realizes. However, its dynamics is quite unusual. The 
point is that although gluons are electrically neutral, their dynamics is strongly influenced by an external magnetic 
field, as one can see from expression (||) for their polarization operator. In a more formal language, while quarks 
decouple and do not contribute into the equations of the renormalization group in that infrared region, their dynamics 
strongly influence the boundary (matching) conditions for those equations at fc ~ m^- 

A conventional way to describe this dynamics is the method of the low energy effective action. By taking into 
account the polarization effects due to the background magnetic field, we arrive at the following quadratic part of the 
low-energy effective action of gluons: 



(2) 

g,eff 



E 

A=l 
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(16) 



where 



a. \e„B\ 



67r ^ 



q=l 1 



1 / \ 1/3 



q—1 ^ y 



> 1. 



(17) 



By making use of the quadratic part of the action as well as the requirement of the gauge invariance, we could easily 
restore the whole low-energy effective action (including self-interactions) as follows: 



r \^ ( f?^ f?^ ^ c TT^A r,A tdA tjA r>A\ 



A=l 



(18) 



where the (chromo-) dielectric constant e = 1 + k was introduced. Also, we introduced the notation for the chromo- 
electric and chromo-magnetic fields as follows: 



Bf = Is,,, {d,At - d,Af + 9f^^^'AfA<^) . 



(19) 
(20) 
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This low energy effective action is relevant for momenta |fc| < nid- Notice the following important feature of the 
action: the coupling g, playing here the role of the "bare" coupling constant related to the scale m^, coincides with 
the value of the vacuum QCD coupling related to the scale i/jei?! (and not to the scale m^). This is because g is 
determined from the matching condition at |fc| nid, the lower border of the intermediate region ^ |fc| <^ ■\/|ei3|, 
where, because of the pseudo-Higgs effect, the running of the coupling is essentially frozen . Th erefore the "bare" 
coupling g indeed coincides with the value of the vacuum QCD coupling related to the scale y^\eB\: g ~ gg. Since this 
value is much less that that of the vacuum QCD coupling related to the scale mj^, this implies that the confinement 
scale Xqcd of the action ( [l8|) should be much less than Aqcd in QCD without a magnetic field. 

Actually, this consideration somewhat simplifies the real situation. Since the LLL quarks couple to the longitudinal 
components of the polarization operator, only the effective coupling connected with longitudinal gluons is frozen. For 
transverse gluons, there should be a logarithmic running of their effective coupling. It is clear, however, that this 
running should be quite different from that in the vacuum QCD. The point is that the time like gluons are now 
massive and their contribution in the running in the intermediate region is severely reduced. On the other hand, 
because of their negative norm, just the time like gluons are the major players in producing the antiscreening running 
in QCD (at least in covariant gauges). Since now they effectively decouple, the running of the effective coupling for the 
transverse gluons should slow down. It is even not inconceivable that the antiscreening running can be transformed 
into a screening one. In any case, one should expect that the value of the transverse coupling related to the matching 
scale TTirf will be also essentially reduced in comparison with that in the vacuum QCD. Since the consideration in this 
section is rather qualitative, we adopt the simplest scenario with the value of the transverse coupling at the matching 
scale rrid also coinciding with g^. 

In order to determine the new confinement scale Aqcd, one should consider the contribution of gluon loops in the 



perturbative loop expansion connected with the anisotropic action (18), a hard problem being outside the scope of 
this paper. Here we will get an estimate of Xqcd, without studying the loop expansion in detail. Let us start from 
calculating the interaction potential between two static quarks in this theory. It reads: 

V{x, y, z) ^ ^\ = (21) 

Because of the dielectric constant, this Coulomb like interaction is anisotropic in space: it is suppressed by a factor 
of ^/e in the transverse directions compared to the interaction in the direction of the magnetic field. The potential 
( pT| ) corresponds to the classical, tree, approximation which is good only in the region of distancies much smaller than 
the confinement radius rgco ~ ^qcd- Deviations from this interaction are described by loop corrections. Let us 
estimate the value of a fine structure constant connected with the perturbative loop expansion. 

First of all, because of the form of the potential ( pT| ) , the effective coupling constants connected with the parallel 
and transverse directions are different: while the former is equal to g^jj.- = gs, the latter is = gs/e^^'^. On the 
other hand, the loop expansion parameter (fine structure constant) is g^j ^ / Ai:Vg, where Vg is the velocity of gluon 
quanta. Now, as one can notice from Eq. (p^), while the velocity of gluons in the parallel direction is equal to the 
velocity of light c = 1, there are gluon quanta with the velocity vj^ = l/-\/e in the transverse directions. This seems to 



suggest that the fine structure coupling may remain the same, or nearly the same, despite the anisotropy: the factor 
in (g^j)^ will be cancelled by the same factor in u^. Therefore the fine structure constant can be estimated as 
as = 5s/47r (although, as follows from Eq. (p^), there are quanta with the velocity = 1, their contribution in the 
perturbative expansion is suppressed by the factor 1 / ^/l) . 

This consideration is of course far from being quantitative. Introducing the magnetic field breaks the Lorentz group 
S'0(3, 1) down to 50(1,1) x 50(2), and it should be somehow manifested in the perturbative expansion. Still, we 
believe, this consideration suggests that the structure of the perturbative expansion in this theory can be qualitatively 
similar to that in the vacuum QCD, modulo the important variation: while in the vacuum QCD as is related to the 
scale |ei?|, it is now related to much smaller scale rn^- 

By making use of this observation, we will approximate the running in the low-energy region by a vacuum-like 
running: 

- ' 6oln-^, where 6o = — (22) 



a'si^i) as 127r 

where the following condition was imposed: a's{md) = as- From this running law, we estimate the new confinement 
scale, 

/a \ 
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We emphasize again that expression (^3|) is just an estimate of the new confinement scale. In particular, both the 
exponent, taken here to be equal to b/bo, and the overall factor in this expression, taken here to be equal 1, should 
be considered as being fixed only up to a factor of order one. 

The hierarchy Xqcd ^ ^qcd is intimately connected with a somewhat puzzling point that the pairing dynamics 
decouples from the confinement dynamics despite it produces quark masses of order Aqc d or less [for a magnetic field 
all the way up to the order of (fO TeV)^]. The point is that these masses are heavy in units of the new confinement 
scale Xqcd and the pairing dynamics is indeed weakly coupled. 



V. QCD WITH LARGE NUMBER OF COLORS 

In this section, we will discuss the dynamics in QCD in a magnetic field when the number of colors is large, in 
particular, we will consider the ('t Hooft) limit Nc oo. Just a look at expression (Q) for the gluon mass is enough 
to recognize that the dynamics in this limit is very different from that considered in the previous sections. Indeed, as 
is well known, the strong coupling constant Ug is proportional to 1/Nc in this limit. More precisely, it rescales as 



(24) 



where the new coupling constant ctg remains finite as Nc — > oo. Then, expression (Q) implies that the gluon mass 
goes to zero in this limit. This in turn implies that the appropriate approximation in this limit is not the improved 
rainbow approximation but the rainbow approximation itself, when both the vertex and the gluon propagator in the 
SD equation are taken to be bare. 

In order to get the expression for the quark in this case, we can use the results of the analysis of the SD equation 
in the rainbow approximation in QED in a magnetic field pO[ | , with the same simple modifications as in Sec. ^ The 
result is 



— C\eqB\ exp 



TvNr 



1/2- 



where the constant C is of order one. As Nc oo, one gets 



C\eqB\ exp 



1/2- 



(25) 



(26) 



It is natural to ask how large Nc should be before the expression (g5|) becomes reliable. From our discussion above, 
it is clear that the rainbow approximation may be reliable only when the gluon mass is small, i.e., it is of the order of 
the quark mass or less. Equating expressions (^ and (p5|), we derive an estimate for the threshold value of Nc- 



N 



thr 



INu + Nd 
In|ei3|/A|p^ 



exp 



2^/l 



II In 



\eB\ 



A 



QCD 



1/2 



(27) 



Expression ( p5| ) for the quark mass is reliable for the values of Nc of the order of N^^^ or larger. Decreasing Nc below 
Nl^^ , one comes to expression (|l|). 

It is quite remarkable that one can get a rather close estimate for iV*'"" by equating expressions (|l|) and (|^): 



thr 



2N^ + Nd 



exp 



11 In. 



\eB\ 



A, 



QCD 



1/2 



(28) 



[notice that the ratio of the exponents ( |27| ) and (|28|) is equal to 0.91]. The similarity of estimates ( p7[ ) and ( |28| ) implies 
that, crossing the threshold iV^*'*'-, expression (Ij) for smoothly transfers into expression (0). 

These estimates show that the value of iV*'"' rapidly grows with the magnetic field. For example, taking Kqcd — 
250 MeV and iV„ = 1, iV^ = 2, we find from Eq. (|^ that A^*'"^ - 10^, 10^, and 10"* for \eB\ - (1 GeV)^, (10 GeV)^, 
and (100 GeV)^, respectively. 



As was shown in Sec. IV, in the regime with the number of colors Nc <^ N^ , the confinement scale Xqcd in 



QCD in a strong magnetic field is essentially smaller than Aqcd- What is the value of Xqcd in the regime with Nc 
being of the order of N^'^^ or larger? It is not difficult to see that Xqcd — Aqcd in this case. Indeed, now the gluon 
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mass and, therefore, the contribution of quarks in the polarization operator are small (the latter is suppressed by the 
factor 1/Nc with respect to the contribution of gluons). As a result, the /3-function in this theory is close to that in 
QCD without a magnetic field, i.e., \qcd — ^qcd- 

Expression (|2^) implies that, for a sufficiently strong magnetic fields, the dynamical mass rriq is much larger than 
the confinement scale Aqcd- Indeed, expressing the magnetic field in terms of the running coupling, one gets 



^QCD exp 



f 



1/2 



(29) 



and for small values of bag ~ A^cQ^s = (i.e., for large values of |ei?|) the mass rrig is indeed much larger than Aqcd- 
This point is important for proving the reliability of the rainbow approximation in this problem. Indeed, the relevant 
region of momenta in this problem is <C \k'^\ <C \eB\ |10| where, because ^ ^qcd ^ strong enough field, 
the running coupling is small. Therefore the rainbow approximation is indeed reliable for sufficiently strong magnetic 
fields in this case. 



VI. CONCLUSION 



QCD in a strong magnetic field yields an example of a rich, sophisticated and (that is very important) controllable 
dynamics. Because of the property of asymptotic freedom, the pairing dynamics, responsible for chiral symmetry 
breaking in a strong magnetic field, is weakly interacting. The key point why this weakly interacting dynamics 
manages to produce spontaneous chiral symmetry breaking is the fact that it is essentially 1+1 dimensional: in the 
plane orthogonal to the external field the motion of charged quarks is restricted to a region of a typical size of the order 
of the magnetic length, £ — 1/ ^y\eqB\. Moreover, such a dynamics almost completely decouples from the dynamics 
of confinement which develops at very low energy scales in the presence of a strong magnetic field. 

While the pairing dynamics decouples from the dynamics of confinement, the latter is strongly modified by the 
polarization effects due to quarks that have an effective 1 + 1 dimensional dynamics in the lowest Landau level. As a 
result, the confinement scale in QCD in a strong magnetic field Xqcd is much less than the confinement scale Aqcd 
in the vacuum QCD. This implies a rich spectrum of light glueballs in this theory. 

This picture changes drastically for QCD with a large number of colors {Nc ^ 100 for \eB\ ^ 1 GeV^). In that case 
a conventional confinement dynamics, with the confinement scale Xqcd ^ Aqcd, is realized. 

The dynamics of chiral symmetry breaking in QCD in a magnetic field has both similarities and important differen- 
cies with respect to the dynamics of color superconductivity in QCD with a large baryon density . Both dynamics 
are essentially 1+1 dimensional. However, while the former is anisotropic [the rotational 5*0(3) symmetry is explicitly 
broken by a magnetic field], the rotational symmetry is preserved in the latter. This fact is in particular connected 
with that while in dense QCD quarks interact both with chromo-electric and chromo-magnetic gluons pO| , in the 
present theory they interact only with the longitudinal components of chromo-electric gluons. This in turn leads to 
very different expressions for the dynamical masses of quarks in these two theories. 

Another important difference is that while the pseudo-Higgs effect takes place in QCD in a magnetic field, the 
genuine Higgs (Meissner-Higgs) effect is realized in color superconducting dense quark matter. Because of the Higgs 
effect, the color interactions connected with broken generators are completely screened in infrared in the case of color 
superconductivity. In particular, in the color-flavor locked phase of dense QCD with three quark flavors, the color 
symmetry is completely broken and, therefore, the infrared dynamics is under control in that case ||2^. As for dense 
QCD with two quark flavors, the color symmetry is only partialy broken down to SU (2)c, and there exists an analog 
of the pseudo-Higgs effect for the electric modes of gluons connected with the unbroken SU(2^c■ As a result, the 
confinement scale of the gluodynamics of the remaining SU{2)c group is much less than Aqcd [ ^7| , like in the present 
case. The essential difference, however, is that, unlike QCD in a magnetic field, the infrared dynamics of a color 
superconductor is isotropic. 

Last but not least, unlike the case of QCD with a nonzero baryon density, there are no principal obstacles for 
examining all these results and predictions in lattice computer simulations of QCD in a magnetic field. 
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